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Motivated by the MIT experiment [Gorlitz et ai, Phys. Rev. Lett. 87, 130402 (2001)], we 
analytically study the effect of density and phase fluctuations on various observables in a quasi one- 
dimensional degenerate Bose gases. Quantizing the Gross-Pitaevskii Hamiltonian and diagonalize 
it in terms of the normal modes associated with the density and phase fluctuations of a quasi-one 
dimensional Bose gas. We calculate dynamic structure factor S{q, w) from low-energy condensate 
density fluctuations and flnd that there are multiple peaks in S{q, to) for a given momentum q due 
to the discrete spectrum. These multiple peaks can be resolved by a two-photon Bragg pulse with a 
long duration which transfer the momentum to the system. We calculate the momentum transferred 
Pz (t) by using the phase-density representation of the Bose order parameter. We also calculate the 
single-particle density matrix, phase coherence length, and momentum distribution by taking care 
of the phase fluctuations upto fourth-order term as well as the density fluctuations. Our studies on 
coherence properties shows that ID Bose gases of MIT experiment do not form a true condensate, 
but it can be obtained by a moderate changes of the current experimental parameters. 

PACS numbers: 03.75.Lm,05.30.Jp 



I. INTRODUCTION 

Recently, there have been a great interest to study 
theoretically [1,2] as well as experimentally [3] on the 
low-dimensional Bose gases. Among the various interest- 
ing topics related to the trapped Bose-Einstein condensa- 
tion (BEG), the study of structure factors and coherence 
properties of quantum Bose gas has attracted major in- 
terest. The dynamic structure factor plays an important 
role which helps us to understand the discrete modes due 
to the finite size system [5]. Unfortunately, the dynamic 
structure factor calculated from the local density approx- 
imation do not show up the discrete nature of the modes 
and this is valid under certain conditions [6,7]. The dy- 
namic structure factor can be analyzed through the be- 
havior of the momentum transferred to a trapped Bose 
gases by two-photon Bragg pulse with long duration [4] 
which reflects the underlying discrete spectrum [5]. On 
the other hand, coherence has importance for many appli- 
cations, like matter wave interferometry, guided atomic 
beams and atom lasers. The realization of the quasi-one 
dimensional Bose systems at MIT [3] raises the question 
of the effects of quantum and thermal fluctuations. The 
discrete nature of the excitations implies that BEG will 
not be destroyed immediately as interactions between 
bosons are turned on [1,2]. In low dimensions, long- wave 
length density and phase fluctuations lead to a new inter- 
mediate state between a condensate and non-condensate 
system. The system forms a quasi- condensate, in which 
the phase of the order parameter is only coherent over 
a flnite distance less than the system size, whereas the 
phase is coherent over a distance of the order of the size 
of the system in a true condensate. Quasi-condensate 
of a phase fluctuating BEGs have been the subject of 
theoretical studies [8-11], including the development of a 



improved many-body T-matrix theory valid in all dimen- 
sions and at low temperatures [12,13], and the extension 
of the Bogoliubov theory to low-dimensional degenerate 
Bose gases [14]. Also, the phase fluctuating nature of 
highly elongated BEG was experimentally demonstrated 
in Ref. [15-17]. 

In this paper, we investigate in details on the prop- 
erties of the effect of the low-energy modes, particu- 
larly, dynamic structure factor, and also the momen- 
tum transferred by a two-photon Bragg pulse since there 
have not been any study in the literature on these prop- 
erties. Moreover, in the theoretical studies on degen- 
erate Bose gases have neglected the contribution from 
the density fluctuations and have concentrated only upto 
the quadratic phase fluctuations. The density fluctua- 
tions and higher order phase fluctuations are important 
to study coherence properties for a wide range of the 
system parameters. In this work, we calculate one-body 
density matrix, phase correlation length and momentum 
distribution by taking care of the density fluctuations and 
fourth-order phase fluctuations. 

This paper is organised as follows. In section II, us- 
ing the phase-density representation of the Bose order 
parameter, we write down an effective hydrodynamic 
Hamiltonian and diagonalize it in terms of the normal 
modes of the density and phase fluctuations. First, we 
rederivc the known low-energy modes at zero tempera- 
ture. Then we incorporate the temperature into theory 
through the temperature dependent chemical potential. 
In section III, we calculate static and dynamic structure 
factors. We also present an analytic expression for the 
time-evolution of momentum transfer to the system due 
to the Bragg pulse. In section IV, we consider the effect 
of the density fluctuations and the fourth-order term in 
the phase fluctuations and calculate the phase coherence 
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function, the phase correlation length and momentum 
distribution. We also discuss the possibility to form a 
true condensate by changing the system parameters. In 
section V, we give a brief summary of our work. 



II. QUANTIZATION OF ID BOSE GAS 

We consider atomic Bose gases confined in a quasi- 
one dimensional harmonic trap as described in MIT lab 
[3]. For all the Figs, in this paper, we will be using the 
parameters used at the MIT experiment on quasi-one di- 
mensional BEC by Gorlitz et al. [3]. In particular, they 
have used ^"^ Na in the trap with axial trapping frequency 
= 27r X 3.5 Hz, {a^ = 1.12 x 10~^ m) and radial trap- 
ping frequency cjq = Stt x 360 Hz, (ao = 1.10 x 10^^ 
m). The three-dimensional s-wavc scattering length is 
a = 2.75 X 10~^ m. The number of condensate atoms 
varies from Nca ^ 10'' to Nco ^ 10'"' which will be speci- 
fied in the Figs. The chemical potential at zero tempera- 
ture is /io = 7iu;^(1.06^7Vco)^/^ ~ 407iu;^ for N^o ~ 10^. 
The parameters satisfy the condition for the system to 
be quasi-one dimensional: TvjJq » /xq >> Tiujz. 

The grand-canonical Hamiltonian of the Bose system 



is 



H' = H-nN 



+ ^ I dz'4}\z,t)i)\z,t)^{z,t)i}{z,t). 



= / dz^\z,t){- — ^ + +-m^lz'' - ii)i){z,t) 



(1) 



The equation of motion for the bosonic field operator 

,^d^{z,t) 



dt 



= [xP{z,t),H'] 

r ^ ^ 2 2 1 7 / \ 

+ gi'^\z,t)'^{z,t)'ii){z,t), 



(2) 



where the effective coupling strength is gx = 2hoJoa. The 
Bose field operator can be decomposed in the as usual 
way, 



^iz,t)^^{z,t) + ij{z,t), 



(3) 



then the equation of motion for the condensate wave 
function is 

.^d^iz,t) . 1 2 2 



+ giMz, t)\H{z, t) + 2<?id$(z, t)h{z, t) 
+ gi^*{z,t)m{z,t). (4) 

Here, n{z,t) = {'ip^z,t)'ip{z,t)) and m{z,t) = 
{ilj{z,t)ip{z,t)) are the normal and anomalous non- 
condensed particle densities, respectively. When tem- 
perature is very low, one can neglect h{z, t) and fh{z, t), 



then the equation reduces to the usual Gross-Pitaevskii 
equation: 



H' = j dz^*{z,t)[- 
I jdz\^{z,t)\\ 



1 



2m5z^ + 2"^"^'^'-^]^(^'*) 
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The condensate density is 



"■0(2:) = — [mo - \rnujlz% 
91 ^ 



(5) 



(6) 



We use the quantized hydrodynamic theory developed 
by Wu and Grifiin [18] at T = 0. The Bose order parame- 
ter can be written as t) = \Jn{z, The fluc- 
tuations in the density and phase about their equilibrium 
aren(.2,f) = nQ{z^-\-^n{z,€) &x\A4>{z,€) = (po{z)+S(f){z,t). 
By keeping up to second-order in the fluctuations, we 
obtain the effective hydrodynamic Hamiltonian at zero 
temperature, 

H = Ho + ^ j dz [mnQ{z)5vl{z, t) + gi5h^{z, t)] , (7) 

where dvz{z,t) = ^j^5(l){z,t) is the superfluid velocity. 
This Hamiltonian is quadratic in the density and the 
phase fluctuation operators, then we can diagonalize it 
using the canonical transformations: 

5h{z, i) = ^ [aji}j{z)e-^i^aj + H.c] , (8) 



and 



5^{z, t) = Y, [hji>i{z)e-"'''aj + H.c] , (9) 



The operators aj and d] destroy and create excitations 
with energy TiLOj and satisfy the commutation relations 
[dj,Q!j,] = [cij,Q;j] = and [Q;j,dt] = 0. Also, 

[5h{z,t),5(j){z' ,t)] — iS{z — z'). The constant terms are 



— i\J ^^'^ = ^ T^- "^^^ equations for the den- 
sity and phase fluctuations can be obtained by using the 
Heisenberg equation of motion: 



d^5n{z,t) 



d_ 
'd'z 



and 



d'^5(j){z,t) d 



—no{z)-^Sn{z,t) 
m oz 



—m{z)^5<i){z,t) 

m oz 



df^ dz 
The equation for the eigenfunctions ^^^{z) is 



(10) 



(11) 
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gid_ 
m dz 



d 

no(.z)—%l>j{z) 



0. 



(12) 



Therefore, Eq.(12) becomes Legendre equation with the 
eigenfrequencies given by ujj = uJz\/3{j + l)/2 and the 
corresponding normaHzed eigenfunctions are i'jiz) = 

\J^^^^Pj{z I Zq) , where Pj(z/Zo) is the Lcgondrc polyno- 
mial in z and Zq is the Thomas- Fermi half length [2] . The 
breathing mode corresponds to the j = 2 which oscillates 
with th(^ fr(^qncncy 0^2 = V^^z ^^nd it has also been veri- 
fied experimentally [19]. The chemical potential depends 
on the number of particles in the condensate and the the 
number of particles in the condensate strongly depends 
on the temperature, Nc{T) = Nco{l - T/Tc), where is 
the critical temperature of an ideal Bose gases [1]. The 
simplest way of including the non-condensate is to use the 
temperature dependent condensate number Nc{T) in the 
chemical potential and therefore in the Thomas-Fcrnii 
half-length. This approximation is similar to the local 
density approximation. Therefore, the chemical poten- 
tial at finite temperature can be written as 



(13) 



and then the Thomas-Fermi half-length at finite temper- 
ature T can be obtain from mulZ^ = 2/i. This tempera- 
ture dependent Thomas- Fermi length, Zt, will be useful 
to study the coherence properties. The condensate nor- 
mal mode frequencies at finite T are the same as at T = 
in the Thomas-Fermi limit, since the frequencies at T = 
do not depend on the value of Nc in this limit, but the 
eigenfunctions changes through the TF half-length. We 
note that the amplitude of the quasi-particle increases as 
temperature increases. 



III. STRUCTURE FACTORS AND BRAGG 
SPECTROSCOPY 

Dynamic Structure Factor. Consider a low-intensity 
off-resonant inelastic light scatters with momentum 
transfer Uq and energy transfer hu to the target (BEC). 
If the external light couples weakly to the number density 
of the target, the differential cross section is proportional 
to the dynamic structure factor S{q,u)), which is ob- 
tained from the Fourier transform of the time-dependent 
density- density correlation functions, 

Siq,Lo) = J J rf^e'('^*-«^) < 5n{z,t)5n{0,0) > . (14) 

It is the density fiuctuation spectrum that can be mea- 
sured in the two-photon Bragg spectroscopy. The dy- 
namic structure factor of this system can be written as. 



X [/,-<5(u; + ujj) + (1 + /,)(5(u; - u;,)], (15) 



where fj = [exp{f3hujj) — is the thermal Bose- 

Einstein function and tpj{q) = J dze~^'^^tpj{z/ZT) is the 
spatial Fourier transform of ipjiz). It can be easily shown 
that 



\^j{qr = n{2j + l)ZT- 



(16) 



where the dimensionless variable is q = qZx- It is conve- 
nient to rewrite the dynamic structure factor as 



Siq,oj) = ^ Sj{q)[fjd{uj + Uj) + {1 + fj)6{u3 - ujj 



where 



sAq) 



4a/2 aaz 



-(2i + i)V7(JTT)- 



(17) 



These functions determine the weight of the light- 
scattering cross-section in S{q,u)) of the corresponding 
collective modes of energy ujj. In Fig.l we have plot- 
ted Sj{q) which as a function of the dimensionless wave 
vector q = qZo for the excitations j — 1,2,3,4, and 5. 
Fig.l shows how many modes significantly contribute to 
S{q,Lu). Note that j = mode do not contribute to 
S{q,oj). It is clear from the Fig.l that the strongest 
weights for these collective modes appear for qZo > 2. As 
an example, for the parameters used at MIT experiment 
[3] which gives Zq ~ 9az, this means that the momen- 
tum transfer g in a light scattering experiments should be 
q > 0.22a^^ for N^o = 10'' in order to pick up the strong 
spectral weight from the low-energy collective modes. We 
also notice that the number of modes that contribute to 
the S{q,oj) increases with the chemical potential. 




FIG. 1. Plots of Sj{q) as a function of the dimensionless 
wave vector q. 
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For q ^ 0, the leading term arises from the Kohn or 
dipole-sloshing mode j = I, with ~ g^; the next con- 
tributions arise from the terms with j = 2, 3. In Fig. 2 
we plot the dynamic structure factor S{q, lo) for the mo- 
mentum transfer corresponds to q = qZ^ = 4. For finite- 
energy resolution we have replaced the delta function by 
the Lorentzian with a width of F = 0.06a;2. 

S(q,w) 




FIG. 2. Plot of the dynamic structure factor S{q, uj) vs. u) 
at T = with momentum transfer q = 4. 

In Fig. 3 we present the results for a higher momentum 
transfer q = 10 using the same arbitrary units as in the 
Fig. 2 , showing multiple peaks at higher energy transfer 
CO. This is expected since for large q, the wave-length 
is small compared to the system size and it excites the 
higher energy modes. 
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FIG. 3. Plot of the dynamic structure factor S{q, uj) vs. u) 
at T = for momentum transfer q = 10. 

As one can see from the Fig. 2 and Fig. 3, the dynamic 
structure factors has multiple peaks and the number of 
peaks increases if the given momentum transferred is 
high. This phenomena is due to the underlying dis- 
crete spectrum. Note that the energy transfer w is much 
smaller than the gap between the chemical potential /i 
and the transverse excited state (~ uq) and hence the 
radial modes can not contribute to the dynamic struc- 



ture factor. This behavior of the multiple peaks can be 
resolved in two-photon Bragg spectroscopy, as shown by 
Steinhauer et al. [4]. 

Bragg Spectroscopy: The observable in the Bragg scat- 
tering experiments is the momentum transferred to the 
condensate. When the system is subjected to a time- 
dependent Bragg pulse, the additional interaction term 
appears in the Hamiltonian which is given by [20], 

Hi{t)= j dzil>\z,t)[VB{t)cos{qz-u)t)]tp{z,t). (18) 

Here, we supposed that the Bragg pulse is switched on at 
time t = and q is also along the z-direction. The mo- 
mentum transferred Pz{t) by the Bragg pulse has been 
calculated based on the Bogoliubov transformation [20]. 
Here, we calculate Pz (t) by using the phase-density repre- 
sentation of the Bosc order parameter. After linearizing 
the Bose order parameter, the above term becomes, 



Hi{t) = J dzno{z)VB{t)cos{qz — ut) + aje 



dztpj{z)VB{t)cos{qz — u)t) + H.c. 



(19) 



The second term in the above equation describe the 
scattering between the condensate and the quasiparticles 
which occurs due to the energy and momentum transfer 
from the optical potential. The time-dependent exponen- 
tials, e^*'^J*, which multiply the quasiparticle operators 
in Eq. account for the free evolution due to the Hamilto- 
nian H' [given in Eq. (5)], and so the Heisenberg equa- 
tion of motion 



becomes 



<') = [{aje-^-^%H + Hj{t)] 



d 



(20) 



(21) 



After solving the equation, 

"jW = Q!j(0) 

dt'V{t')e'^^^' I dzaj^j{z)cos{qz-Lut'). (22) 



The momentum transfer from the optical potential is 

Pz{t) = ^% < a]{t)a,{t) >= (Y^fY^nqSM 



3,g 



X Fj [{^3 - w) , t] - Fj [{ujj +uj),t], 



(23) 



where Fj[{ujj±uj),t] = C^li^'/^JyJ^^ )^ and Sj{q) is given 
in Eq. (17). For large t, S{q,uj) ~ P^(t). We plot P^(i) 
for various time in Fig. 4 and Fig. 5, and shows that the 
multipeak spectrum in S{q, u>) can be resolved only when 
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the duration of the Bragg pulse is t >> 2tt/uJz- When 
t < 2'k/u)z, Pz{t) reflects the dynamic structure factor 
calculated from the local density approximation. 




W/Wz 



FIG. 4. Plots of the momentum transferred Pzit) vs. uo at 
T = when momentum transfer is g = 4 for various time 
U:,t = 7.0 (dashed), uj^t = 12.0 (dotted) and uj^t = 19.0 
(solid). Also, we have assumed Vb = Q-Vhuj^ 
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FIG. 5. Plots of the momentum transferred Pzit) vs. uj 
aX T = Q when momentum transfer g = 10 for various time 
uozt = 7.0 (dashed), u^t = 12.0 (dotted), and uJzt = 19.0 
(solid). Here, Vb ~ O-lhuz- 

Static Structure Factor. For completeness, we also cal- 
culate the static structure factor which can be obtained 
from the Fourier transform of the equal-time density- 
density correlation function 



^Sj{q)coth{ 



iShujj 



(24) 



with the long- wavelength hmit S{q) ^ q^coth{^^) 



IV. DENSITY MATRIX AND MOMENTUM 
DISTRIBUTION 

The single-particle density matrix at equal-time is de- 
fined as 



i?i(z,z') = (^^(z,i)^(z',t)). 



(25) 



The bosonic field operator 'ip{z,t) can be written in 
terms of the phase-density representation as ipiz, t) — 
g»*'(z) ^ h{z). Using the phase-density representation of 
the bosonic field operators into Eq. (25), we get 



Di{z,z') = (v/^e-^(")e*W^')v/^(IO)- 



(26) 



Most of the theoretical studies (except [14]) have con- 
sidered 6n{z) ^ since Sn{z) ^ This is true 
for large interaction strength, but for small interaction 
strength there is a finite contribution to the density ma- 
trix. Therefore, we must include the contribution from 
the density fluctuations to calculate one-body density 
matrix for a wide range of the system parameters. Ex- 
panding square root of the density operator upto second 
order in the following way. 



\/n{z) = \/nn{z) + 6n{z) 



5n{z) 1 5n{z) 2 



-oi: 



2no{z) 8 noiz) 



(27) 



Using the commutation relations for the phase and den- 
sity operators and applying the Wick's theorem, then the 
single-particle density matrix can be written as. 



Di(z, z') = ^no{z)no{z'){{l - 



1 6n{z) 5h{z') 2 



)) 



X e 



8 no{z) no{z') 

i[S4>iz,t)-S(f,{z' ,t')]y ^28) 



We can write down the contributions from the density 
fluctuation to the correlation function as 



((1- 



l,Sh(z) 6n(z'),^,, _ 1 



8 no{z) no(z') 



where 



Sh{z) Sh{z') 
no{z) no{z') 

,5 J" 



^T.i^j + i)VjiITT)cothi^) 



'l-z^ 



1 - z'' 



(29) 



Interestingly, the contribution of the density fluctu- 
ations given in the above is similar to the result ob- 
tained by Mora and Castin [14] by discretizing the low- 
dimensional space. 
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Now we calculate the thermodynamic average of the 
phase part. Following the Ref. [21], we denote, X = 
-i[S^{z,t)-S${z',t')], then 

(e^)=eI^.^l"<^">^ (30) 

The first term and the third term in the exponent is 
zero since (aj) = 0. The second term is, 

The fourth-order term is ^(X*)^ = ^{(^'') - 3(X2)2}. 
We will take the fourth-order term in our following cal- 
culations and show that it contributes significantly to 
the correlation functions even if we use the current ex- 
perimental parameters. Then the single-particle density 
matrix becomes 



Ci (z) 



1 




0.2 0.4 0.6 0.8 1 



FIG. 6. Plots of the normalized one-body density ma- 
trix Ci{z) vs. the the separation z (a) at T = O.STc and 
Fd = Fi^ 0, (b) r = 0.5Tc and Fa = F4, = 0, (c) and (d) at 
T = O.STc and T = O.lTc respectively with a = 2.75 x 10"" 
and other parameters are fixed. 



where 

F.[z,z'] = ^-^ff:^^lP,iz)-P,iz')r 

X cothi^)^ (31) 
and F4[z,z'] = Fi[z,z'] - 3(i^2[^, ^'])^ where 

X (1 + 2/, +2/J). (32) 

We calculate all the summation within the phonon 
regime, hujj < fi and the upper cut-off limit, jmax, is 
obtained from the relation, /i = hojj^^^. The normalized 
one-body density matrix or the phase correlation func- 
tion, Ci(z), is defined as 

^/ nQ(z)no{z') 

This modified one-body density matrix can also be used 
for other quasi-condensates, like quasi-condensate in two- 
dimensional [9] and elongated but three dimensional Bose 
systems [11]. When the system is very close to the 
zero temperature, the contributions from F^i [z\ and [z\ 
are negligible, but these two terms contribute signifi- 
cantly for a large separations and at finite temperature 
(T > 0.2Tc). In Fig. 6, we show the correlation function 
with and without those terms at very high temperature. 



Fig. 6(a) and Fig. 6(b) shows that there is a signifi- 
cant difference between the two curves at large separa- 
tions. To describe the correlation function more accu- 
rately for a wide range of the parameters, we must need 
to include the effect of the density fiuctuations, fourth- 
order term of the phase fluctuations and temperature- 
dependent Thomas-Fermi length Zt- Our studies on co- 
herence properties improves on previous studies in [8]. 
By taking considerations of density fluctuations and the 
phase fluctuations upto fourth-order term, we calculate 
the phase coherence Ci{z) vs. the separation z which is 
shown in the Fig. 7. 



Ci (z) 




z 



0.2 0.4 0.6 0.8 1 

FIG. 7. Plots of the normalized one-body density matrix 
Ci{z) vs. the the separation z for various temperatures. 

The phase coherence length [Lp) of the condensate is 
defined as the distance at which the first order correlation 
hmction is equal to 1/^, i.e., Ci{z — z' = Lp) = 
which is shown in the Fig. 8 for various temperature. 
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0.2 0.4 0.6 0.8 1 

FIG. 8. Plots of the phase coherence length Lp vs. tem- 
perature T in units of Tc- 

The momentum distribution can be obtained from 
Fourier transformation of the one-particle density ma- 
trix: 

n{q)= I ^ dz I ^ dz'Di{z,z')e'''^'-'">. (34) 

Fig. 9 shows that the momentum distribution broadens 
due to the increase of the phase and density fluctuations 
with increasing T. There is a long tail in the momentum 
distribution at high temperature, whereas at low temper- 
ature, it decreases very fast. The long-tail in the momen- 
tum distribution can be used to identify the presence of 
a large phase fluctuations in the bosonic system. 



0.8 - 



0.2 - 








2 4 6 8 10 

qZr 

FIG. 9. Plots of the momentum distribution n(q)/n(0) 
vs. momentum qZr for various temperature T = (dot- 
ted), T = O.lTc (solid), T = 0.3Tc (dashed) and T = 0.5Tc 
(dot-dashed) . 

The estimated temperature of trapped Bose gases at 
MIT is r ~ O.lTc, and Fig7. shows that the phase is 
not coherent overall the system at T = O.lTc. There- 



fore it forms a quasi-condcnsatc with a large phase fluc- 
tuations. One can see from Eqs. (29) and (31) that 
Fd ~ i^m^?'^, and Fa ~ (i^)'/^- It implies that 
for a given number of particles, strong interaction sup- 
presses the density fluctuation where as the phase fluc- 
tuation is enhanced. By moderate change of the inter- 
action strength a, one can produce a true condensate 
even at temperature T ~ 0.1 — 0.27^. In fact, if we use 
a = 2.75 X 10~^^m and keeping other parameters fixed as 
used in other figures, wc find there is a true condensate 
where the phase correlation function Ci (z) is almost con- 
stant (see Fig. 6(c) and Fig. 6(d)) overall the system even 
at temperature T ^ 0.1 — 0.2rc. Also, we plot the mo- 
mentum distribution in Fig. 10 and shows that there is no 
broadening and long-tail in the momentum distribution 
even at T ~ O.BTc. 



0.8 - 
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0.5 1 1.5 2 2.5 3 3.5 

FIG. 10. Plot of the momentum distribution n{q)/n{6) vs. 
qZr at T = O.lTc (dotted) and at T = O.STc (solid) when 
a = 2.75 X 10-"m. 

Note that the new parameters also satisfy the condi- 
tion for a system to be quasi-one dimensional. 



V. SUMMARY 

In this paper, we have presented a detail investigation 
of the effect of density and phase fiuctuations on phys- 
ical observables, like the discrete spectrums, dynamic 
structure factor, momentum transfer due to Bragg scat- 
tering and the momentum distribution. We have also 
studied the dynamic structure factor which reflects the 
discrete nature of the spectrum and discuss how and 
when these multipeaks in the dynamic structure factor 
can be observed through the Bragg spectroscopy. We 
have calculated the one-body density matrix, correlation 
length and the momentum distribution by considering 
the density fluctuations, phase fluctuations upto quartic 
term and also the temperature-dependent Thomas- Fermi 
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length. Our studies improves on the previous studies 
where the density fluctuations, higher order phase fluc- 
tuations and the temperature-dependent Thomas-Fermi 
length were not considered. The coherence properties 
studied from the density and phase fluctuations reveals 
that the Bose gas of MIT experiment do not form a true 
condensate, but we have shown that a true condensate 
can be achieved by a moderate change of the interaction 
strength which are accessible due to the Feshbach reso- 
nance [22]. 



[1] N. J. van Druten and W. Ketterle, Phys. Rev. A 54, 656 

(1996). 

[2] T.-L. Ho and M. Ma, J. of Low Temp. Phys. 115, 61 
(1999). 

[3] A. Gorlitz et al, Phys. Rev. Lett. 87, 130402 (2001). 
[4] J. Steinhauer, N. Katz, R. Ozeri, N. Davidson, C. Tozzo, 

and F. Dalfovo, Phys. Rev. Lett. 90, 060404 (2003). 
[5] C. Tozzo and F. Dalfovo, New J. Phys. 5, 54 (2003). 
[6] F. Zambelli, L. Pitaevskii, D. M. Stampcr-Kurn, and S. 

Stringari, Phys. Rev. A 61, 063608 (2000); A. Brunello, 

F. Dalfovo, L. Pitaevskii, and S. Stringari, Phys. Rev. 

Lett. 85, 4422 (2000); 
[7] D. M. Stamper-Kurn et. al, Phys. Rev. Lett. 83, 2876 



(1999) ; J. M. Vogels et. al, Phys. Rev. Lett. 88, 060402 

(2002) . 

[8] D. S. Potrov, G. V. Shlyapnikov, and J. T. M. Walraven, 

Phys. Rev. Lett. 85, 3745 (2000). 
[9] D. S. Petrov, M. Holtzmann, and G. V. Shlyapnikov, 

Phys. Rev. Lett. 84, 2551 (2000). 
[10] Y. Kagan, V. A. Kasurnikov, A. V. Krasavin, V. N. 

Prokof 'ev, and B. V. Svistunov, Phys. Rev. A 61, 043608 

(2000) . 

[11] D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Wahaven, 

Phys. Rev. Lett. 87, 050404 (2001). 
[12] J. O. Andersen, U. Al Khawaja, and H. T. C. Stoof, Phys. 

Rev. Lett. 88, 070407 (2002). 
[13] U. Al Khawaja et. al, Phys. Rev. A 68, 043603 (2003). 
[14] C. Mora and Y. Castin, Phys. Rev. A 67, 053615 (2003). 
[15] S. Dcttmcr et al, Phys. Rev. Lett. 87, 160406 (2001). 
[16] S. Richard et al, Phys. Rev. Lett. 91, 9053 (2003). 
[17] F. Gerbier, et al, Phys. Rev. A 67, 051602 (2003). 
[18] W-C Wu and A. Griffin, Phys. Rev. A 54, 4204 (1996); 

D. L. Luxat and A. Griffin, Phys. Rev. A 67, 043603 

(2003) . 

[19] H. Moritz et al, cond-mat/0307607 (to appear in Phys. 
Rev. Lett.) 

[20] P. B. Blakic, R. J. Ballagh and C. W. Gardiner, Phys. 

Rev. A 65, 033602 (2002). 
[21] J. W. Kane and L. P. Kadanoff, Phys. Rev. 155, 80 

(1967). 

[22] S. Inouye et al, Nature(London), 392, 151 (1998); J. 
Stenger et al, Phys. Rev. Lett. 82, 2422 (1999). 



8 



